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Abstract: 

The 0{a) corrections to 77 — > ff in the Standard Model are calculated for arbitrary, 
light fermions /. The relevant analytical results are listed in a form that is appropriate for 
practical applications, and numerical results for integrated cross sections are discussed. 
The corresponding QED corrections are generally of the order of some per mille for arbi- 
trary energies. The weak corrections to 77 — > e~e + are negligible below the electroweak 
scale, reach the per-cent level at a few hundred GeV and grow to about —10% at 2TeV. 
The weak corrections to uu and dd production have a shape similar to the one for e~e + , 
but they are larger by factors ~ 1.4 and ~ 3, respectively. 
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1 Introduction 



Since the suggestion of a photon linear collider (PLC) in the 80's [ [3J as an additional 
option for future e + e~ linear colliders, many studies on the feasibility (see Ref. [ |2| and 
references therein) and the physics potential [ 0, [| of such a machine have been performed. 
A PLC provides an excellent device complementary to e + e~ colliders, as can be seen from 
the following examples. Photon-photon collisions allow for a search of Higgs bosons by 
s-channel production and for high-precision tests of the properties of W bosons, which are 
produced in pairs with an enormously large cross section. Moreover, the production cross 
sections of charged particles, in many models for new physics, are even larger than for 
comparable e + e~ machines [ |3j. Last but not least, a PLC allows various QCD studies, 
in particular the investigation of the structure of the photon itself. 

According to the DESY/ECFA study [ a total 77 luminosity of 10 33 cm~ 2 s _1 , or 
even 1-2 orders of magnitude higher, can be reached by Compton backscattering of laser 
photons off the high-energetic beams at a 500 GeV collider. This production mechanism 
renders the luminosity spectrum non-trivial, since both photon beams are not monochro- 
matic, and a luminosity monitor has to be sensitive to both photon energies. For this task 
the processes 77 — > e~e + ,/i~/z + have been suggested (see Ref. [ and references therein) 
as reference reactions. Thus, the lepton-pair production cross section should be known to 
very high precision. 

In this paper we calculate the complete 0(a) corrections to 77 — > ff in the Standard 
Model for arbitrary light fermions /. We present analytical results that are sufficient for 
an evaluation of all relevant observables such as cross sections and distributions for all 
polarization configurations. The structure of the radiative corrections and the leading 
contributions are discussed in detail. Moreover, we provide numerical results on the 
integrated cross sections and the corresponding electroweak corrections for the different 
fermion flavours. 

The paper is organized as follows: in Section we introduce some conventions and 
list analytical results for the lowest-order cross sections. In Section |3] the electroweak 
radiative corrections are classified into QED and weak corrections, and the corresponding 
analytical results are presented. The numerical results are discussed in Section |], and 
Section [5] contains a summary. Explicit analytical expressions for the relevant scalar one- 
loop integrals are given in the appendix. 

2 Conventions and lowest-order cross section 

We consider the reaction 

7 (*!, Ai) + 7(^2, A 2 ) — f(p, a) + j(p, a). (1) 

The mass m/ of the fermion / is neglected whenever possible. Otherwise we follow closely 
the conventions of Ref. [ Q. The helicities of the incoming photons and of the outgoing 
fermions are denoted by Ai ; 2 = ±1 and a, a = ±1/2, respectively. In the centre-of-mass 
system (CMS), the momenta read 

kZ = E(l, 0,0,-1), fc£ = £7(1, 0,0,1), 

p» = £7(1, -sin 0,0,- cos 0), = £7(1, sin 0,0, cos 9), (2) 
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Figure 1: Tree diagrams for 77 — > // 



where E is the energy of the incident photons, and 6 denotes the scattering angle. The 
Mandelstam variables are given by 



s = (h + k 2 ) 2 = (p + p) 2 
t = (k 1 - P ) 2 = (k 2 -p) 2 
(h-p) 2 = (k 2 - P ) 2 



u 



4E 2 } 

-4£ 2 sin 2 §, 
-4£ 2 cos 2 i 



(3) 



The neglect of mj in the kinematics implies that our results are valid for s, —t, —u ^> m 2 . 

The scattering amplitude of 77 — > // obeys Bose symmetry with respect to the 
incoming photons and — neglecting quark mixing — also CP symmetry. Consequently, the 
polarized cross sections d<7 Al,A2 ' <T,cr are related by 



da XlM '^(s, t, u) = d(T A2 ' Al ' <T ' s (s, u, t) (Bose) 
= da- Xi '- x ^- 9 '-°(s,u,t) (CP) 
= do- A2 >- Al -°- a (s,t,u) (Bose + CP). 



(4) 



In lowest order, 77 — >• // is a pure QED process and is therefore invariant under parity. 
Hence, the Born cross sections obey the additional relations 



dff Bom' ,7 ' fT (M,M) 



da 



-Vb — A2, — cr,— a . 
Born 1 



S, t, u) 



(P)- 



(5) 



The two lowest-order Feynman diagrams^ are shown in Fig. [TJ. The differential Born 
cross section reads 



do-Bc 



dQ 



Nf 



647r 2 s , _ 4 

Ai,A2,cr,(T 



£ 7 (l + A 1 P 1 )(l + A a P 2 )|A<^ ff ( a) t ) «)| 2 ) (6) 



where Px,2 are the degrees of beam polarization, and the sum on the r.h.s. includes the 
desired polarizations of the outgoing particles. The colour factor for the fermion / is 



denoted by Nf, i.e. N } 

^Bom ,<T ' CT are g iven b y 



lcpton 



1 and iVq Uark = 3. The squares of the helicity amplitudes 



\M 



Ai ,X2,cr,cr , 
Born ' 



S, t, u) 



AQ\e A - for Ai = -A 2 = ±1, a = -a = ±i 



' t 
t 



4Q 4 /e 4 



u 







for Ai 



otherwise. 



-A, 



=Fl, a 



-a 



2 ■ 



(7) 



1 All Feynman diagrams in this work have been drawn with the help of FeynArts [ 
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The t- and w-channel poles in the squared amplitudes lead to kinematical singularities 
in the very forward and backward directions, where we are not interested in the cross 
sections, since the fermions escape into the beam pipe. For leptons, these singularities are 
of course regulated by a finite lepton mass. For light-quark production in the forward and 
backward directions purely perturbative calculations are not reliable, since the splitting 
of a photon into a nearly collinear quark-ant iquark pair involves QCD effects at very low 
scales. We avoid the forward and backward regions by imposing the angular cut 

cut < e < 180° - e cut . (8) 

For later convenience, we introduce the step function 

9cut(0) = 6(0 - #cut)e(180° - 6 cnt - 9), (9) 

where Q(x) is the usual Heaviside distribution. Integrating over a symmetric angular 
range (|8|), the contributions of all non- vanishing Born cross sections are equal, and the 
integrated, unpolarized cross section reads 



( 1 + COS cut \ 

In — - cos 9 cut 

V 1 - cos 9 cut I 



(10) 



where a = e /(An) is the fine-structure constant. 



3 Electroweak radiative corrections 



3.1 Classification of O(et) corrections and general remarks 

Since 77 — > ff is a pure QED process in lowest order, the SM electroweak corrections 
of 0(a) consist of two separately gauge-invariant types: pure QED corrections and gen- 
uinely weak corrections. The QED corrections include real photon emission (see Fig. |2[), 
virtual photon exchange (see Fig. §), and the corresponding counterterms. The weak cor- 
rections comprise all one-loop diagrams (and contributions to counterterms) that involve 
the massive weak gauge bosons W and Z. For vanishing fermion mass mj, there are no 
contributions involving Higgs-boson exchange or closed fermion loops.0 As a consequence 
the 0(a) corrections do not depend on the Higgs-boson and top-quark masses and on the 
running of a. The weak corrections can be further classified into two subset sQ. The first 
of these subsets includes all diagrams that contain internal Z-boson lines (see Fig. Q), and 
the corresponding corrections are called neutral- current (NC) corrections in the following. 
The second subset includes all diagrams with W-boson exchange (see Fig. f|), leading to 
charged- current (CC) corrections. Note that only the CC corrections involve non-abelian 
couplings among the gauge bosons. 

The perturbative QCD corrections can be obtained from the QED corrections by 
substituting the electromagnetic coupling factor Q"ja by the strong coupling factor 4a s /3. 

2 There are actually Feynman diagrams involving fermion-loop contributions to the AAZ*, AA\* , and 
AAH* vertices. However, these contributions are proportional to the mass of the produced fermion / [ 
H and are thus neglected. 

3 In the i?£ gauges these subsets are gauge- independent. 



3 




+ crossed graphs 

Figure 2: Diagrams for photon bremsstrahlung in 77 // ("crossing" means interchang- 
ing the incoming photon lines.) 




+ crossed graphs 



Figure 3: Diagrams for 77 — > // with virtual Z-boson or photon exchange 




Figure 4: Diagrams for 77 — > // with virtual W-boson exchange 
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The definition of a proper two-jet cross section is, however, problematic, because jets of 
radiated gluons cannot be distinguished from those of quarks. Consequently, a two-jet 
cross section includes the case where a gluon and one of the quarks cause the two jets, 
whereas one of the quarks disappears in the beam pipe. This contribution is divergent, 
owing to the t- or u-channel poles that are not cut out by this definition of two-jet events. 
Therefore, a consistent definition of two-jet events in 77 collisions necessarily involves 
non-perturbative effects, which will not be discussed in this paper. 

The renormalization of the 0(a) corrections turns out to be extremely simple for 
77 — > ff. For rrif = 0, the mass renormalization drops out, and only the wave-function 
renormalization of the external fields and the charge renormalization are relevant. Note 
that the photonic wave-function renormalization constant exactly cancels against a corre- 
sponding part in the charge renormalization (see for instance Ref. [|J) so that no effects 
from the photonic vacuum polarization remain. Consequently, there is no running in the 
electromagnetic coupling a for 77 — > ff in this order. The remaining part of the charge 
renormalization is the contribution of the photon-Z-boson-mixing self-energy at zero- 
momentum transfer, which in the usual 't Hooft-Feynman gauge consists of a W-boson 
loop only and is thus part of the CC corrections. 

Virtual one-loop corrections are included into predictions by replacing the squared 
Born amplitude |A^Bom| 2 by |A^Bom| 2 + 2 Re{A^i_ioo P -A/f Born }, where -Mi_ioo P is the contri- 
bution of the one-loop diagrams to the scattering amplitude. Thus, the one-loop correction 
to a lowest-order cross section is zero whenever the lowest-order vanishes, and we can fac- 
torize the one-loop correction d<7i_i op to the differential cross section into the lowest-order 
cross section do"Born and the relative correction 5i.\ oop for each polarization configuration: 



According to the above decomposition, <5i_i o P is split into NC, CC, and QED contributions. 
Since the Born amplitudes are non- vanishing only for Ai = — A2 and a = —a [see (0)], 
we introduce p = sgn(Ai) = — sgn(A2) and k = sgn(cr) = — sgn(er), and indicate the 
polarization configurations for the relative corrections 5... in flTT|) by S P,K . 

The calculation of the one-loop diagrams has been performed by applying the stan- 
dard techniques summarized in Ref. [ |5| . More precisely, tensor one-loop integrals are 
algebraically reduced to scalar integrals, as described in Ref. [ [7], and scalar integrals are 
computed using the methods and results of Ref. [ . Technically, the algebraic evaluation 
of the Feynman amplitudes, which have been generated with FeynArts [ [5| , has been car- 
ried out in the same way as described in Ref. [ || for 77 — ► tt. In particular, the algebraic 
manipulations have been performed again twice, once using FeynCalc [ |9| and once using 
our own Mathematica [ [l(| routines. For 77 — > e~e + the virtual corrections are related 
to the ones of e~7 — > e~7 [ |TT| by crossing symmetry, which has been used as additional 
check for this channel. 

The evaluation of the real-photonic bremsstrahlung will be described in detail below. 

3.2 Weak corrections 

The NC corrections for the different polarization channels are related by Bose and 
parity transformations as follows: 




^i-ioop — <5nc + <5cc + 5, 



•virt 



(11) 



QED' 
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where gjj- z is the generic Zff coupling, 



9ffz 



s /? 

Cw 5' W C W 



(13) 



with Ij = ±| denoting the weak isospin of the left-handed component of the fermion /. 
The cosine c w of the weak mixing angle is fixed by the ratio of the masses Mw and Mi of 
the weak gauge bosons, i.e. = 1 — = M^/M^. According to (0), all NC correction 
factors can be deduced from 
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a / 



NC 



7T 



+ 



[s + M 7 2 ) 2 



1 



Mf 
u 



'3 u M| 
2 + 1 + ^7 



(t - M 
t 2 



In 



Ml, 



In 



'M z 2 + ^ 



_Mf' 
+ Li 2 



In 1 



u 



Mi, 



Li 2 



m z 2 , 



2\2 

z) 



In — = In 1 - — 7 + Li 2 —9 



Ml 



M z 2 , 



M z 2 , 



Mi_Mi_5\ 
t 2u 4 [ ' v ; 



Note that the contributions to the one-loop correction <5i_i 00 p (|TTJ) are real quantities; the 
imaginary parts of the one-loop integrals do not contribute. 

The CC corrections vanish for right-handed fermions, and the corrections for left- 
handed fermions are related by Bose symmetry, 



%5 = °' 



-'cc 



; cc 



t<->u 



(15) 



For S C q we explicitly obtain 

a 



5, 



cc 



Rei i - 3t " +M >- 2M w t BwW 



+ 



+ 



(Qf-Qf) 2(M 



w 



ut 



u 



Qf 

(Qf-Qf) 2 (2t 



tC ww (t) + uC ww (u) 



Q 2 f 



u 



[B w (t) - B ww {s)) + 



st(t + 2M^)(t-u) 



A, 



[s,t) 



+ 



+ 



t(u-t- 2M 2 V ) 



(M 



w 



u 



u 



WWW 

(s, t) + sC ww (s) + 2tC ww (t) + sC ww (s / 
-st-2t 2 )D www (s,t) 



+ (sM^ - su - 2w 2 )£> www (s, m) - 2sC www (s) 



+ 



QAQr-Qf) 2(M 2 ,- M ) 



Q 2 



f 



(tu + sM^)D ww (u, t) - tC w (t) - uC w (u) 



+ f f| [B(s) _ Bw(t)] _f!(^wz^M CwW 



o 2 , 



+ 



(s + M 2 ^ 2 



s(M^ - t)D w (s, t) + sC(s) + 2tC w (i) 



+ 



(M& - u) 



ir 



s(M^ - u)D w (s, u) + sC(s) + 2mC w 



u 



(16) 



where Q/> = Qf — 21 J denotes the charge of the weak isospin partner /' of the fermion 
/. The functions B , C, .., and D are scalar one-loop integrals, the explicit expressions 
of which are collected in the appendix. Although some of the scalar integrals contain 
(logarithmic) mass singularities, which are regularized by infinitesimal masses rrif and 
ruff, all mass singularities drop out in the final results for 5cc- 

Note that the relative weak corrections vanish directly on the t- and w-channel poles of 
the lowest-order cross sections, i.e. the weak corrections are suppressed where the differ- 
ential cross section is maximal. This is a consequence of the usual charge renormalization, 
which defines the charge e = y/Aira as the 7// coupling for all particles on shell. Since 
this kinematic situation holds for forward and backward scattering, and since the weak 
box diagrams do not develop t- and w-channel poles, all weak corrections are absorbed by 
the corresponding renormalization terms in this special situation. 

3.3 QED corrections 
3.3.1 Virtual corrections 



Bose and P symmetry relate the virtual QED corrections by 



rvirt, p, re 
"QED 



rvirt,— p,K 
"QED 



rvirt,— p, 
"QED 



(17) 



so that it is sufficient to give one particular polarization configuration. Introducing an 
infinitesimal photon mass A <C rrif as IR regulator, and keeping rrif in the mass-singular 
terms, we obtain 



rvirt,+,— 
"QED 



Q 



a 



71 



+ In 



In 



-t 



1 + ln 



m, 



--) + - In [ -- ) + At In 

tj t V u) 2t 2 




The IR divergence drops out after adding the real-photonic bremsstrahlung corrections, 
and the mass singularities cancel completely against mass-singular real corrections caused 
by collinear photon emission, since only final-state particles radiate off photons. 



3.3.2 Real corrections 

Real photon emission in 77 — > ff leads to the kinematically different process 

7 (*!, Ax) + 7(fe, A 2 ) — f(p, a) + f(p, a) + 7 (k', A'), (19) 

with k' and A' denoting the momentum and helicity of the radiated photon, respectively. 
While the incoming momenta k\ t 2 are the same as for 77 — > ff, as specified in @, in the 
CMS, the outgoing momenta read 

= Ef(l, — cos0/ sin9f, — sin0/ sin9f, — cos Of), 
p M = Ef(l,cos(f>f-smQf,sm(f>fSmOf,cosQf), 

k ,,x = cos 0' sinS', sin 0' sin 0', cos 0')- ( 20 ) 
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The lowest-order cross section for 77 — > ffj, which yields an 0(a) correction to 77 — > f f, 
is given by 



N c r 1 

0-7 ( P l> P ^ = ^J dT E 4^ + A l P l)( X + A ^) |>f^^' V | 2 , (21) 



/ dr 

Ai ,A2 ,cr,<7,A 

where the phase-space integral is defined by 



d r = / - f , P / / d ! P / - ^ (2vr) 4 5(A; 1 + k 2 -p-p-k). (22) 
7 (2n) 3 2E f J (2ir) 3 2Ef J (2n) 3 2E' K ' K F y ' y ' 

We have calculated the helicity amplitudes yVl 7 in two different ways. One calcula- 
tion is performed by applying the Weyl-van-der-Waerden spinor technique (see Ref. [ 
and references therein); the second calculation makes use of an explicit representation of 
spinors, polarization vectors, and Dirac matrices. For m/ = the helicity structure forces 
many helicity amplitudes to vanish. In particular, yVf 7 is zero if a = a or Xi = X 2 = —A'. 
Moreover, Bose, CP, P, and crossing symmetry for in- and outgoing photons lead to re- 
lations among the helicity amplitudes. In order to be independent of phase conventions, 
we formulate these relations for |.M 7 | 2 , 



= |/Vl-Ai,-A 2 ,-a,- CT ,-A'|2 £ Cp ^ 



p^p 



J j^-X 1 ,-\2,-<T-a,-X' |2 ^p^j 



I w-A',A 2 ,(r,o-,-Ai 12 



(crossing). (23) 



Owing to these relations, only one independent non- vanishing helicity amplitude is left, 
for which we take 

I .M+'-'-'+'+l 2 = 40 6 ,e 6 (v k i) 2 {p-p) m) 

1 7 1 V/ (p-k 2 )(p-k')(p-k 2 )(p-k'y { ! 

^From this particular |-M 7 | 2 we can read off all different kinds of singularities that can 
occur for 77 — > ff'y. Firstly, there are collinear poles if / or / are scattered into forward 
or backward directions, similar to forward or backward scattering in 77 — > //. In this 
case we again apply angular cuts in order to exclude that / or / escape into the beam 
pipe, i.e. we assume 

9 cut <9 f ,9 f <180° -9 cut . (25) 

Secondly, we encounter the usual soft and collinear singularities if k' becomes soft or 
collinear to p or p. These singularities are the counterparts to the IR and mass singularities 
in the virtual QED corrections given in Section [3.3. It they have to be regularized, as in 
the virtual case, by the infinitesimal photon mass A and the fermion mass nif. In the 
following we describe three different procedures for the treatment of these singularities. 
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(i) IR phase-space slicing and effective collinear factors 

In order to apply phase-space slicing to the IR singularity, we exclude the region 
E' < AE from phase space so that the IR singularity is regularized by the cut energy 
AE <C E. In the soft-photon region A < E' < AE, the asymptotic form of the exact 
differential cross section is known to factorize into the lowest-order cross section without 
photon emission and a universal eikonal factor, which depends on the photon momentum 
(see e.g. Ref. [ §). The integration over the soft-photon phase space, which is carried out 
in the CMS, yields the simple correction factor <5 so f t to the differential Born cross section 
do-Bom for 77 -> //: 



-'soft 




2 In 



2AE\ 



A 




+ iln 2 fa 




(26) 



The factor 5 so f t does not depend on the polarizations of the produced fermions and of 
the incoming photons, and its dependence on A obviously cancels against the one in Sq^ d 
given in Section |3.3.1| . 

The remaining phase-space integration in Q2"T| ) with E' > AE still contains the collinear 
singularities in the regions in which (p ■ k') or (p ■ k') is small. In these regions, however, 
the asymptotic behaviour of the differential cross section (including its dependence on 
nif) has a well-known form (see e.g. Ref. [ |13[). The singular terms are universal and 
factorize from duBorn- A simple approach to include the collinear regions consists in a 
suitable modification of |.M 7 | 2 , which was calculated for mj = 0. More precisely, |-M 7 | 2 
is multiplied by an effective collinear factor that is equal to 1 up to terms of 0(mj/s) 
outside the collinear regions, but replaces the poles in (p ■ k') and (p ■ k') by the correctly 
mass-regularized behaviour. Explicitly, the described substitution reads 



E \m 

A'=±l 



Ai,A 2 ,o-,o-,A' |2 

7 



fAx f ,E f ,a f )f f (x f -,E 



f> a f 



T,f=±l 

x £ \M 

A'=±l 



Al,A2,T<T,T<T,A / |2 

7 I ' 



(27) 



The functions f± describe collinear photon emission with and without spin flip of the 
radiating fermion, 



f + (x f ,E f ,a f ) 
f-{x f ,E f ,a f ) 



' 4£2 s in 2 (f) ' 
4E 2 sin 2 (f )+m? ft 



xA 



Am)E} sin 2 (f ; 



x) + 2x f + 2 [4Ej sin 2 (f ) + m}] 2 ' 



Hit 



where otf = Z(kj,k') is the angle of the photon emission from /. The functions f± 
describing photon emission from / follow by substituting / — > / everywhere. More details 
on this method can be found in Refs. [ |TX|, |I4| , where it is applied to e~7 — > e~77, e'^Z'-f. 

(ii) IR and collinear phase-space slicing 

Instead of using effective collinear factors, one can also apply phase-space slicing to 
the collinear singularities, i.e. the collinear regions are excluded by the angular cuts 



9 



Aa < acf , a,f with Aa <C 1 . The integration over the collinear regions is particularly 
simple for final-state radiation (see also Refs. [[11], Q), since collinear photon emission 
does not affect the kinematics in the factorized Born cross section do"B or n of the non- 
radiative process 77 — > //. The corrections from collinear photon emission can thus be 
described by correction factors S coU to da-Bom, 



da. 



collier, a) 



25+ 11 dcr Bo rn(o-, o) 



c 'colldcrBorn(-0', O ) + C\, oll dcr B orn(o", <J ) 



(29) 



where 



J co\\ 



5. 



coll 




2 In 



(¥) 



5 
2 



2vr 2 



(30) 



Note that the sum of the soft and collinear corrections without spin flip, i.e. <5 so ft + 25+ u , 
comprises all IR- and mass-singular terms originating from real photon emission; after 
adding 5qe D , all In A and mm/ terms drop out. On the other hand, the corrections due 
to S~ on are the only sources for final-state // pairs with a = a. 

(Hi) Subtraction method 

The idea of the subtraction method is to subtract and to add a simple auxiliary 
function from the singular integrand. This auxiliary function has to be chosen such that 
it cancels all singularities of the original integrand so that the phase-space integration of 
the difference can be performed numerically. Moreover, the auxiliary function has to be 
simple enough so that it can be integrated over the singular regions analytically, when 
the subtracted contribution is added again. In the following we apply a modification 
of the so-called "dipole formalism" [ [15] , which was formulated for next-to- leading-order 
QCD corrections involving unpolarized massless partons. In the modified version of this 
formalism all divergences are regularized by photon and fermion masses, and polarization 
is allowed [ [Tj|. 

When the dipole formalism is applied to photon radiation, the combinatorial part in 
the construction of the subtraction function is rather simple. The subtraction function 
consists of contributions labelled by all ordered pairs of charged external particles, one 
of which is called emitter, the other one spectator. Specifically, for 77 — > ff'y we get 
two contributions: in the first case / plays the role of the emitter and / the one of the 
spectator, and vice versa in the second case. The two functions that are subtracted from 
J2\> I -My 1 2 in the phase-space integral are explicitly given by 



\M 



\M 



Xi,\'2,cr,a |2 




2 


sub,l 


(p ■ k') 


l-^i(l 


-yi) 


Ai,A2,o",o-|2 




2 




sub,2 1 — 


(p ■ k') 


.1-22(1 


-2/2) 



-1- Z1 



Z2 



\M 



Ai,A2,c,it, 
Born ' 



S,tl,Ul) 



s,t 2 ,u 2 ) 



(31) 



where lA^g^n' 7 ' 7 ! 2 are the squared Born helicity amplitudes (0) for 77 
variables yi and Zj (i = 1,2) are defined by 

pk' 2pk' 



}Ji 



pp + pk' + pk' 



pp 
pp + pk' ' 



V2 = yi 



ff. The auxiliary 

Z2 = zi _, (32) 



p^p 
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and the Mandelstam variables ti and U{ are defined as in (|3|), but for auxiliary momenta 
Pi and p u 

U = (h - p t ) 2 = (k 2 - p t ) 2 = -AE 2 sin 2 |, 

Ui = (h - Pi) 2 = (k 2 - Pi ) 2 = -AE 2 cos 2 f . (33) 

The auxiliary momenta are chosen such that Pi — > p and pj — > p in the IR limit fe' — > 0, 
that pi — > p + fc' and pi — ► p if k becomes collinear to p, and that p 2 p + k' and 
P2 —> P if k becomes collinear to p. Moreover, the auxiliary momenta obey momentum 
conservation, p + p + k! = pi + pi, and the mass-shell conditions, p 2 = pf = 0, 

Vx - 1 
Pl =p + k - -p, n = p, 

1 - Vi l-yi 

P2 = i P, p 2 =p+k' - 2 p. (34) 

1 — 2/2 1-2/2 

^From this definition we can also deduce that the scattering angles 8i, which are defined 
in (|33|), are given by 9\ = 9j and 9 2 = 6f. 

It is straightforward to check that J2i |-M S ub,;| 2 has the same asymptotic structure as 
J2x' l-M 7 | 2 m t ne s °ft limit k! — > and in the collinear limits (p • k'), (p ■ /c') — > 0, so that 
the phase-space integral 



* = -5f/dr 



7 1 ) 9cut 



(35) 



^ /dr 



7 1 ) fi'cut (9f)9cut(8f) - ( l-^sub^P^cutC^/) ) - ( |-M sub , 2 | 2 fl'cut(6 ) /) 



2s 

is finite and can be performed numerically. In (|35"D we indicated the phase-space cuts 
(p5|) explicitly by the step functions (||). In the subtracted part the cuts are applied 
to the auxiliary momenta pi, pi. Since these approach the physical momenta of the 
final-state fermions in the singular regions, the cuts do not obstruct the cancellation of 
the singularities in (^) as long as they avoid the singularities. The last equality holds 
because 9 2 = 9f and Q\ = 6j in our case. 

For the full cross section we have to add the integral of J2i \-M S ub,i\ 2 that is evaluated 
with the regulators A and rrif [ [RJ. The functions |-M S ub,i| 2 are constructed such that 
the integration over the photon phase space can be performed analytically, leading to 
universal correction factors 5^ b on the Born cross section (XBorn ( |1 0|) for 'yy — > ff, 

a 2 {a,a) = 2<5+ b Cr B orn(o", a) + <^ b £7Born(— C, &) + 5 s 7 lb (TBorn(0", -O"), (36) 

where 

e ub = Q 2 ^<|ln(-]ln(^]+ln 

^sub 






In ( p6[) the Born cross sections 0"B O rn are evaluated with the restriction (|| 
tering angle 9: 8 cut < 9 < 180° — 8 cnt - As required, the IR- and mass-singular terms in 
25^ nh exactly cancel against those terms in 5qe D . The final result for the real-photonic 
contribution to the cross section is given by <r 7 = <j\ + o 2 . 
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4 Numerical results 

For the numerical evaluation we have adopted the parameters [ |l7j 

a = 1/137.0359895, M w = 80.41 GeV, M z = 91.187 GeV. (38) 

We need not specify the masses m/ of the light fermions, since these are only kept as 
regulating parameters and drop out in all considered observables. We discuss only un- 
polarized cross sections. The non-vanishing cross sections for polarized initial states and 
unpolarized final states differ from the unpolarized cross sections only by the normaliza- 
tion. 

In Fig. H we show the lowest-order and the (9(a)-corrected cross section for 77 — > e~e + 
for the angular cuts 8 cnt = 5°, 10°, 20°, 40°. The Born cross sections vary from 137 pb to 
33 pb for these cuts at ^/s = 100 GeV; they scale like 1/s if the cut angle 9 cut is chosen 
energy-independent, as can be seen in (|TI]). Since the impact of the 0(a) corrections is 
hardly visible in Fig. |5], we show the relative QED and weak corrections to 77 — > e~e + 
separately in Fig. |G| for two angular cuts. 

For an energy- independent angular cut 8 cnt , the QED corrections (see Fig. |^) do not 
depend on the scattering energy for s 3> m e 2 , since all electron-mass singularities cancel, 
and s is the only scale that survives. The cancellation of all potentially large QED 
corrections such as a\n(m^/ ' s)/n implies that the resulting QED correction is of the order 
of 0(a/7t), i.e. of the order of several per mille. The numerical results confirm this 
expectation. The weak corrections stay below 0.05% for energies below 100 GeV and tend 
to zero in the low-energy limit. In other words, weak-boson exchange decouples below 
the electroweak scale. Above 100 GeV, the weak corrections become sizeable and develop 
a peak at ^Js = 2M\v, originating from diagrams with a W-pair cut in the s channel. For 
energies up to 1 TeV, £ wea k reaches several per cent and becomes more and more negative 
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77 — > uu, dd 
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y/sf GeV 


$cut 


4m/ P b 


<%ed7% 


^weak/% 


^wcak/% 


^weak/% 


10 


5° 


13722 


1.30 


0.00 


0.00 


0.00 


10° 


10130 


0.74 


0.00 


0.00 


0.00 


20° 


6595.2 


0.33 


0.00 


0.00 


0.00 


40° 


3270.9 


-0.02 


0.00 


0.00 


0.00 


100 


5° 


137.22 


1.30 


0.02 


0.05 


0.26 


10° 


101.30 


0.74 


0.02 


0.07 


0.34 


20° 


65.952 


0.33 


0.03 


0.10 


0.49 


40° 


32.709 


-0.02 


0.05 


0.14 


0.73 


500 


5° 


5.4889 


1.30 


-0.97 


-1.39 


-3.03 


10° 


4.0520 


0.74 


-1.29 


-1.85 


-4.08 


20° 


2.6381 


0.33 


-1.78 


-2.62 


-5.97 


40° 


1.3084 


-0.02 


-2.47 


-3.79 


-9.23 


1000 


5° 


1.3722 


1.30 


-2.81 


-3.88 


-7.95 


10° 


1.0130 


0.74 


-3.61 


-5.03 


-10.49 


20° 


0.65952 


0.33 


-4.70 


-6.69 


-14.48 


40° 


0.32709 


-0.02 


-5.95 


-8.78 


-20.11 


2000 


5° 


0.34306 


1.30 


-6.18 


-8.27 


-15.92 


10° 


0.25325 


0.74 


-7.62 


-10.35 


-20.38 


20° 


0.16488 


0.33 


-9.33 


-12.98 


-26.60 


40° 


0.081773 


-0.02 


-11.15 


-15.98 


-34.55 



Table 1: Integrated Born cross section for 77 — > e e + , the corresponding relative QED 
and weak corrections, and the weak corrections to 77 — > uu, dd 

with increasing energy, crossing the —10% mark at about 2TeV. The large negative 
corrections at high energies are due to Sudakov-type logarithms such as aln 2 (M^/s)/7r, 
and the dominant contributions stem from W-boson exchange. 

The lowest-order cross sections and the relative QED corrections for 77 — > // with 
arbitrary fermion flavour can be easily obtained from the results on 77 — > e~e + by multi- 
plying the results for the e~e + pair by the factors NjQj and Q 2 j, respectively. In particular, 
this means that the QED corrections to non-leptonic channels become even smaller. The 
weak corrections, however, depend on the fermion flavour in a non-trivial way. Therefore, 
the relative weak corrections are explicitly shown in Fig. |6] also for up-type and down- 
type light quarks. The shape of the weak corrections to light-quark-pair production is 
qualitatively similar to the one for lepton-pair production, but they are larger in size. For 
high energies we roughly get 5™^/5q E 6 d ~ 1.4 and ^weak/^QED ~ 3. This enhancement 
of the relative weak corrections is mainly due to the suppression of the lowest-order cross 
section by the quark charges, which is not present in the dominating CC corrections. 
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Method 


AE/E 


Aa/rad 


e C nt = io° 


#cut = 20° 


IR slicing and 
effective collinear factor 


1(T 3 


- 


0.798 ±0.016 


0.345 ±0.014 


KT 5 


0.819 ±0.029 


0.329 ±0.024 


IR and collinear slicing 


1(T 3 


ID" 3 


0.756 ±0.011 


0.3302 ±0.0083 


KT 5 


0.784 ±0.015 


0.349 ±0.013 


KT 5 


io~ 3 


0.734 ±0.019 


0.323 ±0.015 


io~ 5 


0.808 ±0.027 


0.324 ±0.022 


Subtraction scheme 






0.74447 ±0.00080 


0.33124 ±0.00069 



Table 2: Comparison of results for the QED correction 5qed/% at s/s = 500 GeV, ob- 
tained by the different methods for bremsstrahlung corrections described in Section ISTB 



Table [I] summarizes the discussed results by providing some representative numbers. 
At high energies the weak corrections are dominated by the CC corrections. The NC 
corrections are at the level of 1% at 2TeV. 

We conclude our numerical discussion by a short comparison of the different methods 
for the singular phase-space integration for real photon emission, which are described 
in Section |3.3.2| . Table |2] compares numerical results on <5qed that have been obtained 
by performing the multidimensional integration with Vegas [ using the same Vegas 
parameters for each integration. The subtraction method leads to an integration error 
that is smaller by a factor of 10-20 with respect to the results from the two versions of 
phase-space slicing. While there are still large compensations between the phase-space 
integral and the (semi-) analytically calculated singular parts in the slicing approach, for 
the subtraction method all compensations take place between 5qe D and 2<5+ b , which 
are computed without delicate numerical integrations. Table ^| illustrates the consistent 
application of the different methods, but the numbers for the smaller cut # cut = 10° also 
reveal that Monte Carlo integration by Vegas tends to underestimate integration errors 
if the integrand becomes complicated, although it has been smoothed by appropriate 
transformations of the integration variables^. This distinguishes the subtraction method 
that is less sensitive to numerical uncertainties. 



5 Summary 

The 0(a) corrections to 77 — * // in the Standard Model have been calculated for 
arbitrary light fermions /, i.e. fermion-mass effects are neglected. Compact analytical 
results for the cross sections have been listed for arbitrary polarization configurations, 
rendering their incorporation in computer codes very simple. Numerical results on the 
corrections to integrated cross sections have been discussed. 

The corrections are classified into QED and purely weak corrections. Owing to the 
cancellation of all mass-singular contributions between virtual and real-photonic correc- 



4 Repeated evaluations for 6> cut = 10° show that the results obtained with the two slicing variants come 
closer and closer to that of the subtraction method given in Table @ if the statistics is improved. 



15 



tions, the QED corrections to integrated cross sections are of 0(Q"ja/7T) for all energies, 
i.e. of the order of some per mille. For lepton-pair production the weak corrections are 
negligible below the weak bosons scale, reach some per cent at 1 TeV, and reduce the cross 
section more and more with increasing energy, crossing —10% at about 2 TeV. For up- 
and down-type quarks the weak corrections to the integrated cross sections show the same 
qualitative features as in the leptonic case, but the corrections are a few times larger. The 
weak corrections vanish whenever the differential cross sections develop t- or u-channel 
poles, i.e. the relative corrections can be enhanced or suppressed by appropriate angular 
cuts. 

The smallness and the structure of the corrections to 77 — > e~e + ,/i _ /z + underline the 
suitability of these processes as a luminosity monitor. In particular, the corrections do not 
exhibit large uncertainties due to hadronic effects in the photonic vacuum polarization 
or due to the less precisely known top-quark or even unknown Higgs-boson mass. The 
results for the processes 77 — > qq provide a valuable input for QCD studies. 



Appendix 

List of scalar integrals 

Here we list all scalar one-loop integrals that are needed for the evaluation of the 
virtual corrections given in Section [| We use the same definition of the momentum-space 
integrals and of the arguments of the standard functions B , Co, and D Q as given in the 



appendix of Ref. [ |TTJ. The relevant integrals are calculated for the limit |s|, \t\, \u\, ^> 



mj,mj,. By definition, Mandelstam variables with a hat get an infinitesimal imaginary 
part ie, with e > 0, i.e. s = s + ie etc. After supplying this imaginary part where necessary, 



all scalar integrals can also be obtained from those for Compton scattering in Ref. [ |Tl|] 
upon using crossing symmetry. Scalar functions that are related by the interchange of t 
and u are given generically with the abbreviation r = t,u. 

All needed 2-point functions B are calculated in D space-time dimensions with D —>■ 4. 
Instead of using B directly, we have preferred to introduce the UV-finite combinations 



So (s, 0,0)- B (Q, 0,M, 



w 




5o(r,0,Mw)-S (0,0,Mw) = 

B Q (s, M w , M w ) - B Q (0, 0, M w ) = B ww (s) = (3 w \n(x w ) + 1, (A.l) 
with the abbreviations 



The relevant 3- and 4-point functions are given by 
Co(0,0,s,mfi,mfi,mfi) = C(s) — — In 
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C (m 2 f ,0,r,Mw,mf,mf>) = C w (r) 



Co(0,0,s,0,M w ,0) = C w (s) 



Li 2 



Ml 



Li, 1 



-In 



M w . 



M w -r 



7T 

IT 



In 1 - 



M w , 



C (0,0,r,0,M w ,M 



\\j - C w (r) - --Li 2 f jj^r) , 



'w, 



C (0, 0, s, M w , 0, M w ) = C w (s) = - ln 2 (x ¥ 



Co(0,0,s,M w ,Mw,Mw) = C www (s) = — ln 2 (x w ) 



(A.3) 



D o (0, 0, 0, 0, s, r, mf,M-w,mfi,mfi) = D w (s,: 



s(r - M w ) 



Li 2 1 + 



M w , 



4Li 2 



-In 2 

2 V mj, 



s * + 21n [-- i-llnfl 



r \ 7r 



Ml) 6 



-Do(0, 0, 0, 0, t, u, Mw, Mw, m//, m/') = D ww (t,u) 



tu- M w (w + t) 



2 Li 2 ( 1 + x 4u - -^-1^ ] + 2?7 ( -x tu , 1 - ^^r I In I 1 + x tu - ttt-' 



M w 



2Li 2 (l + x 4u ) + ln 



'M w -A ln / _ t 



m), 



M w , 



with x fu 



M w (t + w) 



tM-M w (t + w) 
^(0, 0, 0, 0, s, r, M w , 0, M w , M w ) = D www (s, r) 
1 



E(-i)" +1 



s 2 (f -M w ) 2 -4f 2 sM w „=i 
3Li 2 (l+x„)-Li 2 + 



M 



'M w -f, 



In 1 + 



M w -r / 



+ 111 1 ~ 772~ M -JC n) - J2 ( Li 2 (! + X nXl) + T) {-X n , X*) In (1 + X n X*)} 



M w 



=±1 



(A.4) 



with x w and /3 W as given in (|A.2| ) and 



xi, 2 



s(r - M w ) - 2fM w ± t/s 2 (r - M w ) 2 - 4f 2 sM| 



w 



/[2(f + s)M w ]. (A.5) 
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The dilogarithm Li 2 (x) and the function T](x,y) are defined as usual: 

dt 



f x dt 

Li 2 (x) = - / — ln(l - t), -n < arc(l - x) < 7r, (A. 6) 

Jo t 

rj(x, y) = hi(xy) — ln(x) — ln(y), — n < arc(x), arc(y) < n. (A. 7) 
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